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Abstract 



O , We study the gravitational behaviour of a spherically symmetric radiating star when 

^ ■ the fluid particles are in geodesic motion. We transform the governing equation into 

a simpler form which allows for a general analytic treatment. We find that Bernoulli, 
Riccati and confluent hypergeometric equations are possible. These admit solutions 
in terms of elementary functions and special functions. Particular models contain the 
£^ ■ Minkowski spacetime and the Friedmann dust spacetime as limiting cases. Our infinite 

CN ' family of solutions contains specific models found previously. For a particular metric 

^ I we briefly investigate the physical features, derive the temperature profiles and plot 

^ . the behaviour of the casual and acasual temperatures. 

> 

X 1 Introduction 

Relativistic models of radiating stars are useful in the investigation of cosmic censor- 
ship hypothesis, gravitational collapse with dissipation, formation of superdense mat- 
ter, dynamical stability of radiating matter and temperature profiles in the context 
of irreversible thermodynamics. The general model, incorporating all necessary phys- 
ical requirements and variables, is complicated and difficult to solve; the treatments 
of Herrera et al [1] and Di Prisco et al [2] involving physically meaningful charged 
spherically symmetric collapse with shear and dissipation illustrate the complexity of 
the processes. To solve the field equations, and to find tractable forms for the gravi- 
tational and matter variables, we need to make simplifying assumptions. De Oliviera 
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et al [3] proposed a radiating model in which an initial static configuration leads to 
collapse. This approach may be adapted to describe the end state of collapse as shown 
by Govender et al [4]. In a recent treatment Herrera et al [5] proposed a model in 
which the form of Weyl tensor was highlighted when studying radiative collapse with 
an approximate solution. Maharaj and Govender [6], Herrera et al [7] and Misthry et 
al [8] showed that it is possible to solve the field equations and boundary conditions 
exactly in this scenario. For recent treatments involving collapse with equations of 
state and formation of black holes see Goswami and Joshi [9], [10]. 

A useful approach in understanding the effects of dissipation is due to Kolassis 
et al [11] in which the fiuid trajectories are assumed to be geodesic. In the limit, 
in the absence of heat fiow, the interior Priedmann dust solution was regained. This 
solution formed the basis for many investigations involving the physical behaviour 
such as the rate of collapse, surface luminosity and temperature profiles. These include 
the analytic model of radiating spherical gravitational collapse with neutrino fiux by 
Grammenos and Kolassis [12], the model describing realistic astrophysical processes 
with heat fiow by Tomimura and Nunes [13], and models undergoing collapse with 
heat fiow as a possible mechanism for gamma-ray bursts by Zhe et al [14]. Herrera 
et al [15] considered geodesic fiuid spheres in coordinates which are not comoving 
but with anisotropic pressures. Govender et al [16] showed that the behaviour of the 
temperature in casual thermodynamics for geodesic motion produces higher central 
temperatures than the Eckart theory. The first exact solution with shear, satisfying 
the boundary conditions, was obtained by Naidu et al [17] by considering geodesic 
fiuid trajectories. Later Rajah and Maharaj [18] extended this treatment and obtained 
classes of models which are nonsingular at the centre. 

It is clear that the assumption of geodesic motion is physically acceptable and 
has been used by other investigators in attempts to describe realistic astrophysical 
processes. In this paper we attempt to perform a systematic treatment on the governing 
equation at the boundary for shear-free collapse by assuming the geodesic motion of 
the fiuid particle. Our intention is to show that the nonlinear boundary condition 
may be analysed systematically to produce an infinite family of exact solutions. In 
Section 2, we present the model governing the description of a radiating star using 
the Einstein field equations together with the junction conditions. We show that it 
is possible to transform the junction condition to a Bernoulli equation and a Riccati 
equation. Solutions are obtained in terms of elementary functions in Section 3. In 
Section 4, we show that the boundary condition, under relevant assumptions, can be 
written in the form of a confiuent hyper geometric equation. We demonstrate that an 
infinite family of solutions in terms of elementary functions are possible. In Section 5, 
we obtain the explicit form for the causal temperature using the truncated form of the 
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Maxwell-Cattaneo heat transport equation for a particular metric. This illustrates that 
the simple forms for the gravitational potentials obtained in this paper are physically 
plausible. Some concluding statements are made in Section 6. 



2 The model 

We analyse a spherically symmetric relativistic radiating star undergoing shear-free 
gravitational collapse. This assumption is reasonable when modelling a radiating star 
in relativistic astrophysics. If we suppose that the particle trajectories are geodesic 
then the acceleration vanishes. Then the line element, for the matter distribution 
interior to the boundary of the radiating star, is given by 

ds^ = -dt" + [dr^ + r2(rf^2 ^ gi^2 Q^^2y^ 

where B = B{r, t) is the only surviving metric function. The energy momentum tensor 
including radiation for the interior spacetime is given by 

Tab = {p + P)UaUb + P9ab + QaUb + QbUa (2) 

where the energy density p, the pressure p and the heat flow vector q are measured 
relative to the timelike fluid 4-velocity u"" = 6q. The heat flow vector takes the form 
q"' = (0, q, 0, 0) since q ■ u = for heat flow which is radially directed. 

The nonzero components of Einstein field equations, for the line element ([1]) and 
the energy momentum tensor ([2]), can be written as 

B B^ 1 fB'"^ 2B'\ 
' = -'B-B-^^B-^[-B^^-r^)^ ^2^) 

B B^ 1 fB" B'^ 1B'\ 
" = -^B-W^ + W^y-B-B^+rB)^ ^^'^ 

2 I B' B'b\ , 

where dots and primes denote differentiation with respect to time t and r respectively. 
Equating (l3b|) and (15c|) we obtain the condition 

which is the condition of pressure isotropy. Equation (jlj) is integrable and we obtain 
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where Ci{t) and C2{t) are functions of time, and dis a constant. As the functional form 
for the potential B is specified the matter variables p,p and q are known quantities, 
and the system (3) has been solved in principle. 

The interior spacetime ([1]) has to be matched across the boundary r = b to the 
exterior Vaidya spacetime 



ds' 



2m(f ) 
R 



dv^ - 2dvdR + R^idO"^ + sin^ edcf)), 



(6) 



where m{v) denotes the mass of the star as measured by an observer at infinity. The 
hj^ersurface at the boundary is denoted by S. The matching of the line elements ([1]) 
and ([6]), and matching of the extrinsic curvature at the surface of the star, leads to a 
set of equations. The boundary conditions at S have the form 



dt 



2m dRV^^ , 



Rj:, 

3 



B'B 



B^ 
B 



-r^B' 



(7a) 

(7b) 
(7c) 

(7d) 



where the subscript means that the relevant quantities are evaluated on S. 
From ([3]), ([5]) and fITcl) we generate the condition 



Adh{CiC2 - CiC2)iCib^ - C2) - ACiC2iCib^ - C2f 
-2d\C,b^ - C2)iC,b^ - C2) + 5d\C,b^ - C2Y 



0. 



Effectively results from the nonvanishing of the pressure gradient across the hyper- 
surface S. Equation ([8]) governs the dynamical evolution of shear-free radiating stars 
in which fluid trajectories are geodesic. To complete the description in this particular 
radiating model we need to explicitly solve the differential equation ([H]). 



3 Generating Analytic Solutions 

A particular solution to ([8]) was found by Kolassis et al (1988) by inspection. We 
show that it is possible to transform ([H]) into familiar differential equations which 
admit solutions in closed form. Our method is a more systematic approach in solving 
equation ([8]). In this approach we let 

C^b"" -C2 = u{t). (9) 
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On substituting ([9]) into ([8]) we can write 

Abdu'^Ci + A{u^ - hdu)uCi - Ah'^u^Cl = d^{2uu - 5u^). (10) 

Equation (10) is simpler than ([8]) and can be viewed as a first order differential equation 
in the variable Ci. In general, ffTOj) is a Riccati equation (in Ci), and is difficult solve in 
the above form without simplifying assumptions. For the integration of fllOj) . in terms 
of elementary functions, we consider the following two cases: 



3.1 Bernoulli equation 

We set 



2uu — 5u^ = 
so that the function u is given by 

u = a or u = /?(t + 7)-2/^ (12) 

where a,P and 7 are real constants. With the assumption ffTTl) . ffTOj) becomes 

Abdu^Ci + 4{u^ - bdu)uCi - Ab'^u^Cl = 0. (13) 

Equation f|T3!) is nonlinear but is a Bernoulli equation which can be linearised in general. 
When u = a, equation (fT3l) becomes 

b 



a 



(14) 



which is a Bernoulli equation with constant coefficients. The solution of ([T^ is given 
by 



a 



6^ — exp ^ 



a{t+e) \ ' 



bd 



where e is the constant of integration. Consequently the remaining function C2 has 
the form 



a exp 



a(t+e) 
bd 



Hence the interior line element ([T]) has the specific form 

^2 52 _ 



ds^ = -dr + 



ajt+e) 
bd 



— exp 



f ait+e) \ 



[dr^ + r\de^ + sin^ ed(f)'')] 



(15) 



in terms of exponential functions. We believe that this is a new solution to the Einstein 
field equations for a radiating star. It is interesting to observe that if we set a = c? 
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when t — s> oo (or large values of the constant e) then (flSjl becomes the flat Minkowski 
spacetime 

ds^ = -df + dr^ + r\de^ + sin^ edcj)"^) 

which is a limiting case. 

When u = /5(t + 7)-2/3, ([TSD becomes 



(16) 



which is also a Bernoulli equation with variable coefficients. The solution of ( fT6l) is 
given by 



-2/3 



where / is the constant of integration. Consequently the remaining function C2 is given 
by 

bd 



Co 



( 3/3(t+7)i/a 
I hd 



-2/3 



62 + /3/exp 

Hence the interior line element ([T]) takes the particular form 



ds^ = -dr + 



If we set 



then (fT71) becomes 



62 + /5/exp 



3/3(t+7) 
bd 



1/3 \ n 



+ /3/exp 



/ 3/3(t-H7)i/3 
\^ fed 



sm 



(17) 



1/3 



7 = 0,rf= — bJ = —,P = -- 



6 y 



62 



1/3 



1 - ab^ exp (§) 
l-ar2exp(f)'/' 



n 2 



t4/3 [c^r2 + r2(d^2^sin2e(i(/)2)] 



which was first found by Kolassis et al [11]. Here we have shown that their model 
found by inspection arises naturally as a solution of a Bernoulli equation. It is easy to 
see that for large values of the constant / we obtain 

ds" = -dt" + \^dr^ + ^2^^^2 ^ gi^2 ^^^2)] 

from ( |T7|) . This corresponds to the Friedmann metric when the fluid is in the form of 
dust with vanishing heat flux. 
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3.2 Riccati equation 

If we set 

then the function u is given by 



u — dhii = 



u = —bd(t + a 



,-1 



where a is a constant. In this case equation (|T0|) becomes 

AbdCi - Ah'^Cl + d\t + a) -2 = 0, 



(19) 



(20) 



which is an inhomogeneous Riccati equation. The solution of equation fl20|) has the 
form 



-d 



1 - V2 + (1 + V2)g{t + a)^ 



{t + a)-\ 



(21) 



2b [l+g{t + a)^] 

where g is the constant of integration. Consequently the remaining function has the 
form 

l-V2 + {l + V2)g{t + a)'^ 



C2 = hd{l- 



2[l+^(t + a)v^] 
Hence the interior metric ([T]) has the specific form 

d'^it + af 



ds' 



-dr + 



[Ci{r^ - b^){t + a) - bdf 



[dr^ + r^de^ + sin" 



(22) 



which is written in terms of Ci. We believe that fl22l) is a new solution for a radiating 
star whose particles are constrained to travel on geodesies. The simple form of fl22l) 
will assist in studying the physical features of our model. The solution (12^ arises 
in a natural way once we realise that the underlying dynamical equation at the 
boundary is a Riccati equation. 



4 Special functions 

The solutions found in the previous sections all have power law forms for the quantity u. 
It is possible that other solutions in terms of elementary functions or special functions 
may exist with a power law representation for u. Consequently in this section we 
attempt to generate a general class of solutions to the model dH]) by assuming 

M = a(t + a)". (23) 
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On substituting (l23l) into (ITOl) we obtain 



Lao 



n 



d 



{t + a)Ci - -it + a)^C^ = -7Tn(3n + 2 



(24) 



The nonlinear equation ( l24l) is a Riccati equation and it is difficult to solve the equation 
in the above form. If we introduce a transformation 



h U 

then ([23]) becomes the second order linear differential equation 



{t + afU 



a 
db 



{t + a 



— n 



{t + a)U -U = 



(25) 



(26) 



in the function U with variable coefficients. We can transform fl26l) to simpler form if 
we let 

,„,i i, (n + l)± JAnin + 1) + 1 

Then (|26|) becomes 



2(n+ 1) 



(27) 



(n + l)z^' 



If we let 



#2 

X 



+ 



z/' + 2fc(n + l) 



hd 
—aip 



dW ^^yy g 

dip hd 



(28) 



6(i(n + 1) 
then (l28l) has the equivalent form 



F(X) = Vr(^) 



d'^Y dY 



A;y = 0. 



(29) 



Observe that (!29|) is the confluent hypergeometric equation with solution in terms of 
special functions in general. 

Note that the solution of (!29l) can be written in terms of 

Jik,2k;X), 



Y 



W 



J k,2k; 



-aip 



hd{n + l] 

where J are Kummer functions. In general the solution of the equation fl28|) can be 
written in terms of the Kummer series. Observe that when > we can write 



W 



J{k, 2k- X) 
T{2k) 
fik)f 



e^^ \t(1 



T)f''dT 



(30) 



as a particular solution of the differential equation (!28l) where T{z) = e^'^r^^^dr is 
the gamma function. From fl30l) we note that the solution can be expressed in terms 
of elementary functions for all natural numbers k. Consequently the differential equa- 
tion (12^ admits solutions in terms of elementary functions when k is a natural number. 
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4.1 Particular metrics 



We can regain previous cases from the general form (1501) . We illustrate this feature for 
particular values of k. When k = 1, we obtain n = or n = —2/3. For this case the 
particular solution of the equation (!28l) becomes 



W 



e^-1 



X 



X 



bd{n + 1) 



(31) 



with the help of (130 



When n = 0, from (l27j) and fl3T|) we can easily see that 



U 



bd 



a 



1 — exp 



—a(t + a) 
bd 



is a particular solution of the equation (!26l) . Then with the help of (l25ll we find that 



62 



(32) 



is a particular solution of ( IMl) which is given by 

+ - ^Ci^ = 0. 
bd d 

The general solution of (!33|) becomes 



(33) 



52 



D + exp 



a{t+a) 
bd 



where D is an arbitrary constant. Consequently the interior metric ([T]) has the specific 
form 

-■ 2 

b'd' 



ds^ = -dr + 



D + exp 



bd 



JJ^2 _|_ ^2 gjjp ^ 



Q(t+a) 
bd~ 



[dr^ + r^{de^ + sm^ed4)'')\ (34) 



in terms of exponential functions. Note that the line element (IMI) reduces to the metric 
(fT5l) if we set D = -fe^. 

When n = —2/3, from (l27j) and (I3T1) we observe that 



3a 



1 — exp 



-3a(t + a 



|l/3 



is a particular solution of the equation fl2^ . Hence with the help of fl2^ we obtain 

a(t + a)2/3 



62 



1 — exp 



3a(f+a)i/3 
bd 



(35) 



as a particular solution of (p4|) which has the form 

Ci + 

The general solution of ( l36l) becomes 



Ci - -Ci^ = 0. 



(36) 



aD{t + a' 



2/3 



62 



D + exp 



3a(t+a)V3 
bd 



where D is an arbitrary constant. Consequently the interior metric ([T]) takes the 
particular form 



-1 2 



Dr'^ + exp 



3Q(t+a)V^ 
bd 



sm 



(37) 



Note that the line element (1371) reduces to the line element (fTTl) if we set D 



a/- 



4.2 A new solution 

It is possible to generate an infinite family of new solutions from the general form (1301) 
by specifying values for the parameter k. These may correspond to new solutions for a 
radiating sphere which are not accelerating. We illustrate this process by taking k = 2 
(so that n = —2 or n = —4/5) in ( l30l) . We consider only the case n = —2 as the 
integration procedure is same for other values of k (or n) . For this case the particular 
solution of the equation (128|) becomes 
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W^= — [2 + X + (X-2)e^], X 



When n = —2, from ( l27j) and ( |38l) we observe that 



bd{n + 1) 



(3^ 



U 



6b^d^ 



a'' 



[2bd{t + a) + a] - [2bd{t + a) - a] exp 



a 



bd(t + a) 



is a particular solution of the equation (l26l) . Hence with the help of fl25l) we obtain 

(39) 



2b'^d^(t + a)2 - [2bd(t + a)(M(t + a) - a) + a^] exp 



bd{t+a) 



62 



{2bd{t + a) — a) exp 



fe(i(t+a) 

is a particular solution of (12^ which has the form 



{2bd{t + a)+a) {t + a 



2d 



^it + a)-2 + 2{t + a)-i - -Ci^ = - — (t + a)-2. 
od \ d b 



(40) 
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The general solution of (HOj) becomes 



[a^ + 2bd{t + a){bd{t + a) - a)] exp [^jf^] +2lfdfD{t + a 



62 



D{2hd{t + a) + a) + {2hd{t + a) - a) exp 



(t + a) 



(41) 



M(t+a) 

where D is an arbitrary constant. Consequently the interior metric ([1]) has the specific 
form 

rf2 



ds' 



-dr + 



[dr^ + r^de^ + sin" ed(l)^)] , (42) 



[Ci{r^ - b^) + a{t + a)-^f 

where Ci is given by ( 14T|) . Hence we have found a new solution to the boundary 
condition ([8]) by specifying a particular value for the parameter k. This process can be 
repeated for other values of k and an infinite family of solutions are possible in which 
the gravitational potentials can be expressed in terms of elementary functions. 



5 Physical analysis 



The simple forms of the gravitational potentials found in this paper permit a detailed 
study of the physical features of a radiating star. In this study we consider the partic- 
ular line element flM|) and set a = bd and a = to obtain 

2 



ds' 



-de + 6^ 



D + exp(t) 



[dr^ + r\de^ + sin' 



Dr"^ _l_ ip. exp(t) 

for simplicity. For the metric fHSl) the matter variables can be written as 
3D exp(t) {exp(t) [Db'^ + Qb'^Dr'^ + Dr^ + 46^ exp(t)] + AD\^} 



P 



P 



[D + exp{t)f [Dr^ + b"^ expit)f 



Dexp(t) 



[D + exp(t)]^ [Dr^ + b"^ exp{t)f ' 
{exp(t)) [2&2exp(t)(r2 - 36') - 2Db 

4:Drexp{t) 



2 2 

r 



(43) 
(44a) 

30} , 
(44b) 

(44c) 



[D + exp(t)]" 

When D = then (143!) becomes the Minkowski metric with p = p = q = The 
matter variables are expressed in simple analytic forms which facilitate the analysis of 
the physical behaviour. From (jHj) we have that at the centre of the sphere 

3D [D + 4exp(t)] 



Po 



[D + exp{t)f 
?,D[D + 2exp{t)] 
[D + exp(t)]2 
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so that po and po have finite values at the centre r = with vanishing heat fiux go- The 
gravitational potentials in fH5]) are finite at the centre and nonsingular in the stellar 
interior. The quantities p,p and q are well behaved and regular in the interior of the 
sphere, at least in regions close to the centre. At later times as t ^ oo we note that 
g oc r so that the magnitude of the heat fiux depends linearly on the radial coordinate. 

Next we briefiy consider the relativistic effect of casual temperature of this model. 
The Maxwell-Cattaneo heat transport equation, in the absence of rotation and viscous 
stresses, is given by 

rK'qb + qa = -f^ {h^.T + Tua) , (45) 

where hah = dat + UaUb projects into the comoving rest space, T is the local equilibrium 
temperature, k (> 0) is the thermal conductivity and r (> 0) is the relaxation time. 
Equation (H5ll reduces to the acausal Fourier heat transport equation when r = 0. For 
the line element ([1]), the casual transport equation (145!) can be written as 

T{t,r) = -^J T{qB)B + qB^ dr (46) 

for geodesic motion. Martinez [19], Govender et al [16] and Di Frisco et al [20] have 
demonstrated that the relaxation time r on the thermal evolution, plays a significant 
role in the latter stages of collapse. For the line element fH3l) . fj46|) becomes 

T(t,r) : 



t¥ exp(t) {2D'r' - ¥ exp(t) [exp(t) - D]} 
K [exp(t) + D] [62 exp(t) + Dr'^f 
262exp(t) 



K [6^ exp(t) + Dr"^ 



hit), (47) 



where h{t) is a function of integration. For simplicity we assumed that r and k are 
constant. The function h{t) may be related to the central temperature Tc{t) by 

From flTTI) and fHHl) the temperature can be written as 

rDr"^ {Dr^ [D - exp(t)] - 2b'^ exp{t)} 



T(t,r) = TJt)- „ 

k[D + exp(t)] [Dr^ + exp(t)]' 



K [Dr^ + 62 exp(t)] ' 



(49) 



When r = 0, we can regain the acausal (Eckart) temperature profiles from fH^ . In 
Fig. 1, we plot the casual (solid line) and acasual (dashed line) temperatures against 
the radial coordinate on the interval < r < 5 for particular parameter values [k = 
T = l,b = 5, D = 70 and h{t) = 0) on the spacelike hypersurface t = 1 . We observe 
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that the temperature is monotonically decreasing from centre to the boundary in both 
casual and acasual cases. It is clear that the casual temperature is greater than the 
acasual temperature throughout the stellar interior. At the boundary S we have 

^s)casual ~ '"s)acasual" 

Our figures have been generated by assuming constant values for the parameters r 
and K. Changing the values of the relaxation time and the thermal conductivity would 
produce different gradients for the curves but the result would not change qualitatively. 




Causal 



Acausal 



Figure 1: Temperature T vs radial coordinate r (r = 1). 



6 Discussion 

It is possible to introduce shear in geodesic motion as shown by Naidu et al [17] and 
Rajah and Maharaj [18] in the description of a radiating star. The solutions that follow 
are governed by a Riccati equation and have a complicated form. Consequently in this 
paper we have considered the simpler case of a shear-free metric with particles traveling 
on geodesic trajectories. The master equation, governing the boundary condition of 
the stellar model, was transformed to a simpler form. Under certain assumptions a 
Bernoulli equation is possible. This Bernoulli equation admits two solutions in terms of 
elementary functions: the first solution contains the Minkowski spacetime as a limiting 
case and the second solution corresponds to the Kolassis et al [11] model with the 
Friedmann dust spacetime as the limiting case. A general class of solutions are possible 
if we transform the master equation to a confiuent hypergeometric equation. The 
resulting transformed equation admits solution in terms of special functions namely 
the Kummer functions. By specifying particular values for a parameter in the special 
function we demonstrate that an infinite family of solutions, in terms of elementary 
functions, are possible. The simple form of the solutions makes it possible to study the 
physical features of the model and to find an analytic form for the causal temperature. 
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